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ABSTRACT 

We have observed the spectrum of Procyon A (F5IV) from 4559 to 5780 A with a S/N of ~ 10 3 and 
a resolving power of 2 x 10 5 . We have measured the line bisectors and relative line shifts of a large 
number of Fe I and Fe II lines, comparing them to those found in the Sun. A three-dimensional (3D) 
hydrodynamical model atmosphere has been computed and is tested against observations. The model 
reproduces in detail most of the features observed, although we identify some room for improvement. At 
all levels, the comparison of the 3D time-dependent calculations with the observed spectral lines shows 
a much better agreement than for classical homogeneous models, making it possible to refine previous 
estimates of the iron abundance, the projected rotational velocity, the limb-darkening, and the systemic 
velocity of the Procyon binary system. The difference between the iron abundance determined with the 
3D model and its ID counterpart is < 0.05 dex. We find consistency between the iron abundance derived 
from Fe I and Fe II lines, suggesting that departures from LTE in the formation of the studied lines 
are relatively small. The scatter in the iron abundance determined from different lines still exceeds the 
expectations from the uncertainties in the atomic data, pointing out that one or more components in 
the modeling can be refined further. 

Subject headings: convection — line: formation — line: profiles — stars: abundances — stars: 
atmospheres 



1. INTRODUCTION 

Widely available grids of model atmospheres are based 
on a classical set of assumptions that have remained largely 
unchanged for many years. In brief, the classical assump- 
tions are that the atmosphere consists of plane-parallel 
homogeneous layers in hydrostatic and local thermody- 
namic equilibrium through which a constant (radiative 
plus convective) flux is transmitted. This is not to say 
that model atmospheres have not been progressively re- 
fined; one notes, for example, the increasing completeness 
with which line blanketing is included as an opacity source. 
The geometric assumption of plane parallel homogeneous 
layers has been replaced by the adoption of spherically ho- 
mogeneous layers for construction of model atmospheres of 
supergiants and giants. More recently, the assumption of 
local thermodynamic equilibrium (LTE) has been replaced 
to a limited extent by that of non-LTE, largely for models 
built to represent individual stars. Another classical as- 



sumption - hydrostatic equilibrium - has more tenaciously 
resisted refinement. This paper discusses models that dis- 
pense with this assumption (and the geometric assumption 
of homogeneous layers) by simulating the surface convec- 
tion as a time-dependent hydrodynamic phenomenon. 

The models representing atmospheric convection are a 
development of three-dimensional hydrodynamical mod- 
els that realistically represent solar granulation (Stein & 
Nordlund 1998; Asplund et al. 2000a). Previous ex- 
tensions of the models to stars have been discussed by 
Dravins & Nordlund (1990) (Procyon A, a Cen A and B, 
and (3 Hyi) and Asplund et al. (1999) (metal-poor stars 
HD 84937 and HD 140283). In this paper, we confront the 
models with an extensive set of line profiles of Procyon A 
obtained at very high-resolution and high signal-to-noise 
ratio. Our tests are similar to those discussed by Dravins 
& Nordlund but provide a much stiffer challenge of the 
models. 

Crudely, the 3D hydrodynamical model atmosphere 



1 Departamento de Astrofi'sica, Universidad de La Laguna, E-38206 La Laguna, Spain 

1 



2 



comprises areas of hot rising gas (granules) and cool sink- 
ing gas (inter- granular lanes). Relative to a line profile 
from a classical atmosphere, the lines are broader, shifted, 
asymmetric, and of a different equivalent width. Deriva- 
tion of fundamental parameters, abundances, and other 
quantities extracted from observed spectra has not been 
attempted using 3D stellar model atmospheres with a few 
exceptions (e.g. Bruls & Rutten f992; Atroshchenko & 
Gadun f994; Kiselman & Nordlund 1995; Shchukina & 
Trujillo Bueno 2001). This may be a serious omission, 
because modeling the effect of the convective motions 
and temperature inhomogeneities on the equivalent widths 
through the parameters known as micro- and macrotur- 
bulence may represent an important source of systematic 
errors. To the extent that equivalent widths are used to 
derive defining fundamental parameters, classical and 3D 
hydrodynamical models will yield different results. For 
the solar case, however, Allende Prieto et al. (2001) have 
shown that the abundances derived using a numerical in- 
version based on a ID model atmosphere are accurate to 
within 0.04 dex. Naturally, this will not hold in general 
for other stars and elements, and it is unclear whether the 
same result is valid for classical, theoretical, ID models. 

Observational scrutiny of surface convection has con- 
centrated on the measurement of the asymmetry of lines, 
which is generally quantified as a line bisector (Kulander 
& Jefferies 1966). Line bisectors have been analyzed in 
the Sun (e.g. Dravins, Lindegren & Nordlund 1981), and 
stars across the Hertzsprung-Russell diagram (e.g. Gray 
1982; Gray & Toner 1986; Gray & Nagel 1989). The mea- 
surement of line shifts has not received the same attention 
(e.g. Allende Prieto & Garcia Lopez 1998; Hamilton & 
Lester 1999). We refer the reader to Dravins (1999) for a 
recent review on line asymmetries and shifts generated by 
surface convection. 

Our goal is to test a new 3D model of Procyon using 
a large collection of lines observed at high-resolution and 
high signal-to-noise ratio. Tests employ the line shifts, 
asymmetries, and widths, or in other words, the detailed 
line profiles. This the first such comprehensive test of 3D 
hydrodynamical models for a star other than the Sun. All 
the calculations described in this paper assume LTE. 

2. OBSERVATIONS 

The optical spectrum of Procyon was observed from Mc- 
Donald Observatory on January 30 and 31 (UT), 1999. 
We made use of the Harlan J. Smith 2.7m telescope and 
the 2dcoude spectrograph (Tull et al. 1995), at the fo- 
cal station Fl, delivering a FWHM resolving power of 
R = X/8X ~ 2 x 10 5 . Five different spectral setups, se- 
lected by tilting the E2 grating, provided almost complete 
spectral coverage from 4559 to 5780 A, and different expo- 
sures (up to ten for a single setup) were coadded to reach 
signal-to-noise ratios per pixel in the range 550 — 2000. 

A large number of flatfield exposures was acquired in 
order to correct the pixel-to-pixel sensitivity pattern (typ- 
ically smaller than 0.5 %) without degrading the poten- 
tially high signal-to-noise ratio of the spectra. The indi- 

2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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vidual exposures of Procyon were shorter than, or equal 
to, five minutes to keep the signal within the linear range 
of the detector, and to minimize the blurring due to varia- 
tions in the Earth's motion and instrumental drifts. Stan- 
dard procedures (within IRAF 2 ) were used to remove the 
bias level and scattered light, to correct the pixel-to-pixel 
variations of sensitivity, and to extract the spectra. The 
spectral coverage is complete, except for ten gaps span- 
ning from 0.2 to 12.5 A wide. For each spectral setup, 
about a hundred Th-Ar emission lines were used to define 
the (dry standard air) wavelength scale, by fitting a 2D 
polynomial with an r.m.s. at least ten times smaller than 
the dispersion. The flux distribution along any given order 
was normalized by a polynomial. While the normalization 
procedure is expected to be adequate in spectral regions 
where the density of absorption lines is low, and the lines 
present are not very strong, the limited coverage of the 
individual spectral orders (~ 20 A) precludes an accurate 
continuum tracing in the proximity of strong lines and in 
crowded spectral regions. 

The multiorder spectra permit a precise correction of the 
velocity variations and instrumental shifts between differ- 
ent exposures, cross-correlating them order-by-order. In 
all cases the standard deviation of the shifts measured be- 
tween two frames from the 16 different orders was less than 
50 m s^ 1 , with a mean value of 18 m s _1 (SEM 3 : 4.5 m 
s" 1 ). For a given spectral setup, we shifted the individual 
extracted spectra with a single velocity to overlap them 
on top of the first exposure. Then, the five different spec- 
tral setups were reduced to a common scale by comput- 
ing and correcting the heliocentric velocity shifts produced 
by Earth's rotation, and the motion of the Earth-Moon 
barycenter around the Sun. The extracted orders were 
then combined to compose a spectroscopic atlas, which 
spans the aforementioned spectral range at a variable dis- 
persion of 0.0083 < AA < 0.0109 A/pix. Figure |l| displays 
the spectrum of Procyon in the vicinity of the Mg I line at 
5528.4 A, and compares it with the solar spectrum (Ku- 
rucz, Furenlid & Brault 1984) and the photographic atlas 
of Procyon by Griffin & Griffin (1979). The new data are 
publicly available through the internet 4 . 

In a work of this nature, it is very important to check for 
possible instrumental effects that might distort the stellar 
line profiles. The 2dcoude spectrograph is a well-tested 
instrument but instrumental drifts have been suspected 
(see, e.g., Allende Prieto et al. 1999a). Figure [2] shows the 
variation of the heliocentric correction for the individual 
exposures respect to the first (filled circles), and the ve- 
locity drifts determined from cross-correlating the individ- 
ual exposures with the first as template. A few relatively 
large drifts (up to 50 m s _1 ) are present between consec- 
utive exposures, but slower drifts are also obvious in some 
cases. Spectra of the Th-Ar hollow cathode were acquired 
immediately before and after each stellar exposure. The 
reference wavelengths for the Th-Ar lines used in the cali- 
bration correspond to dry standard air (760 mmHg and 15 
°C), and the profusion of calibration spectra should pre- 
vent large drifts as pressure, temperature, and humidity 




drag the air refraction index. Nonetheless, the sensitiv- 
ity to changes in temperature, pressure and humidity is 
such that variations on small time scales - between the 
exposures of the Th-Ar hollow cathode and the star - 
can definitely be a source of small errors. For example, 
a change of 0.1 °C from standard conditions would induce 
a change of 28 m s _1 , a change of 0.1 % in pressure would 
induce a drift of 82 m s^ 1 , and introducing 0.1 % of wa- 
ter vapor would produce a change of 13 m s _1 at 5000 
A. Lacking a proper study of these effects in our spectra, 
we note that changes in pressure within the spectrograph 
follow in full amplitude those outside the building, and 
variations as large as than 0.07 % per hour were mea- 
sured on the observing dates. Additionally, it is known 
that the position of the CCD varies with the weight of 
the N2 dewar. Differences between how the slit is illumi- 
nated by the Th-Ar hollow cathode and the stellar light 
cannot be ruled out as contributor to systematic errors. 
Although the cross-correlation technique makes it possi- 
ble to determine the shifts among the individual exposures 
for a given setup to within ~ 5 m s _1 , the observed sys- 
tematic effects introduce an additional uncertainty in our 
wavelength calibration, probably not larger than 20-30 m 
s _1 . Considering more than 200 lines in the Th-Ar spec- 
tra, we confirmed that the nominal resolving power was 
indeed achieved (R = 197700 ± 500; see Fig. |), with no 
detectable variation with order or wavelength, but a small 
dependence on the setup. 

3. OBSERVED SIGNATURES OF CONVECTION 

Different factors may contribute to the observed line 
asymmetries and shifts. Blends with other features affect 



particular lines in different ways that can not be gener- 
ally predicted with precision. Therefore, it is necessary to 
average observed line asymmetries and shifts for a large 
number of lines in order to quantify the signatures of con- 
vection with minimal distortion from other shifts. Con- 
vective line asymmetries and shifts are closely related and 
carry complementary information. However, the limited 
spectral range of earlier studies prevented the study of 
line shifts in depth. In short, stellar studies have ignored 
line shifts and relied on mean line bisectors computed as- 
suming that the velocity of the lowest end of the bisector 
is the same for all lines. 

Previous studies of line asymmetries in Procyon's spec- 
trum were limited in scope. Gray (1981a, 1982) analyzed 
lines in a very small spectral window (~ 70 A) acquired 
at R = 1.2 x 10 5 . Rice & Wehlau (1984) attempted the 
observation but at insufficient spectral resolution. Dravins 
(1987) analyzed 11 lines using the ESO coude spectrom- 
eter double-pass photoelectric scanner at R = 2 x 10 5 , 
and he strengthened his results by demonstrating that, 
despite limited signal-to-noise ratio, the scanned version 
of the photographic atlas of Griffin & Griffin (1979) could 
be used to derive statistically-meaningful averaged bisec- 
tors for sets of lines with different strengths. Our much 
larger spectral coverage than previous CCD studies, higher 
signal-to-noise ratio and resolution, and thorough control 
of possible systematic effects, set the stage for the first de- 
tailed study of asymmetries and shifts for a large sample of 
iron lines. The gravitational redshift and the orbital and 
radial velocities affect the zero point of the velocity scale, 
but the information we are searching for is contained in 
the line bisectors and their relative displacements. 
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Fig. 2. — Velocity shifts determined from the cross-correlation of the individual exposures for each setup (solid line with error bars). The 
error bars are determined as the standard error of the mean from the scatter for all the spectral orders. The velocity shifts produced by 
Earth's motion are marked with filled circles. All shifts are relative to the first spectrum in each series. 
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Fig. 3. — Number of emission lines measured in the Th-Ar hollow cathode spectra as a function of the resolving power R = X/SX. Cosmic 
rays are easily identified and excluded from the analysis. Systematic errors from blended features are minimized by fitting only the right flank 
of the nearly Gaussian distribution (thick part of the histogram). 



To minimize systematic differences in the comparison 
with the Sun, we have carried out an identical analysis 
with the solar atlas of Kurucz et al. (1984). The lines were 
selected (or rejected) and classified following the same cri- 
teria for both stars, in an attempt to keep the comparison 
as meaningful as possible. It is of great interest to ana- 
lyze the asymmetries of Fe I lines in the flux spectrum of 
the Sun on an absolute scale, as this has never been done 
consistently before, and must provide the best standard 
for comparisons with other stars. Finally, even at risk of 
repeating what has already been said many times in the 
literature, we note that there are several good reasons to 
pick Fe I among the observed species: the large number 
of lines in late-type stars, the lack of measurable isotopic 
shifts and hyperfine splitting, and the availability of accu- 
rate laboratory wavelengths. 

3.1. Line shifts 

We used the line list of Thevcnin (1989, 1990) to search 
for Fe I lines in the solar spectrum. The central wave- 
lengths and bisectors of the Fe I lines in that list were 
measured in the solar atlas and the McDonald spectrum 
of Procyon, and then compared with the precise laboratory 
wavelengths measured by Nave et al. (1994). The wave- 
length calibration of the solar atlas has been examined by 
Allende Prieto & Garcia Lopez (1998), who showed that 
it is highly reliable. 

The measurement of the central wavelength of the line 
profiles was carried out fitting a third-order polynomial to 
the 11 lowest points around the minimum, equivalent to 
an interval of 90 — 120 mA in the spectrum of Procyon 
and 50 mA in the solar spectrum. We measured all the 
lines in Thevenin's list identified as Fe I lines. For this 
comparison, the equivalent widths of the lines were esti- 
mated by synthesizing the lines with Kurucz (1992) model 
atmospheres corresponding to the Sun and Procyon. This 
is safer than attempting to measure the equivalent widths 
in the spectra, as a large fraction of the lines are strongly 
blended outside their core, especially at the bluer wave- 



lengths. We made use of the solar gf— values of Thevenin 
(1989, 1990) and the damping enhancement factors for Fe 
I recommended by the Blackwell group, as implemented 
in the LTE synthesis program MOOG (Sneden 1973). 

The left panels in Fig. show the Fe I line velocity shifts 
(v) measured in the spectra of the Sun and Procyon. The 
smaller wavelength coverage available for Procyon (4559 
to 5780 A) in comparison with the solar atlas of Kurucz 
et al. (2960 to 13,000 A), reduces the number of mea- 
sured lines from 1551 to 506. The line shifts exhibit sim- 
ilar patterns for both stars. There is no reason to think 
that the well-known granulation observed in the Sun is not 
present in Procyon, and we naturally associate the larger 
shifts towards the red for stronger lines with a decreasing 
of the convective blue-shifts. Unfortunately, there is no 
line in the observed part of Procyon's spectrum stronger 
than 200 mA, leaving open the question of whether the 
velocity shifts of the strongest Fe I lines fall on a solar-like 
'plateau'. It is interesting to compare the left-hand pan- 
els in Fig. ||. The velocity shifts in the solar spectrum 
span a range of ~ 0.6 km s _1 , but at least 0.9 km s _1 for 
Procyon's spectrum. 

The velocity indicated by the strong lines on the solar 
plateau (W\ > 200 mA) is essentially independent of the 
line's equivalent width. We note that Fig. |^ has been 
truncated at W\ — 250 mA, but the solar plateau extends 
at least to lines as strong as 2 A (see Allende Prieto & 
Garcia Lopez 1998). Part of the observed scatter should 
be attributable to the laboratory wavelengths and uncer- 
tainties in the measurement of the center of a spectral line 
(with a finite width), which can be quantified. Nave et 
al. (1994) classify their wavelengths into four categories, 
depending on their uncertainty (en), which ranges from 
0.4 mA to more than 10 mA. Most lines have errors be- 
low 5 mA and about half of the lines have errors below 
1 mA. Errors in the measured stellar wavelengths can be 
estimated from the horizontal scatter of the polynomial 
least-squares fit to the line center divided by the square 
root of the number of points in the fit (in this case 11), o~ s . 
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Fig. 4. — Left panels: velocity shifts between the stellar and laboratory wavelengths for relatively weak lines; for Procyon the zero-point 
arbitrary. The strong solid line is a least-squares linear fit for all lines with 20 < W\ < 180 mA. Right panels: difference between the measured 
shifts and the linear fits to lines with 20 < W\ < 180 mA (histogram); the expected distribution, based on error estimates for the laboratory 
wavelengths and the polynomial fits to the center of the stellar lines, is also shown (solid curve). 



Neglecting the effect of unidentified line blends and 
systematic effects in the wavelength scales of the stel- 
lar spectra, the uncertainties in the measured line shifts, 
a = y/af + af, or more precisely, the distribution of un- 
certainties, n(a), can be used to predict the distribution 
of errors (S) in the observed velocity shifts as 




where C is chosen such that J^ D oo N(S)dS = 1. The de- 
rived N(S) can be compared to the observations in order 
to constrain the intrinsic scatter. 

Discarding a deviant line, the average shift for the 
19 measured Fe I solar lines stronger than 400 mA is 
v = 0.626±0.011 (SEM) km s -1 , consistent with the gravi- 
tational redshift of solar photospheric light at Earth, 0.633 
km s _1 . The strengthening of the convective blueshifts for 
weaker lines is definitely not far from linear. A linear re- 
gression for lines with equivalent widths 20 < W\ < 180 
mA provides, in the solar case, a similar slope as linear fits 
within several subintervals included in that range. Con- 
versely, there is some marginal indication that the trend 
flattens for the lines with 20 < W\ < 60 mA in the spec- 
trum of Procyon. The intrinsic scatter for the Sun is very 
small; at least significantly smaller than our measurement 
errors (see Fig. For Procyon, there is an additional 
scatter, but it should be kept in mind that we have ne- 
glected other sources of systematic errors, such as those in 
the theoretical equivalent widths. We will return to this 
issue in the light of 3D model atmospheres. 



Making use of the laboratory wavelengths measured re- 
cently by S. Johansson (private communication), we derive 
the velocity shifts for 25 Fe 11 lines, which exhibit a trend 
similar to that for Fe I lines. The velocity shifts of Fe II 
lines in the spectrum of the Sun are generally more shifted 
to the blue than the Fe I lines, in agreement with the re- 
sults obtained by Dravins, Larsson, & Nordlund (1986) for 
the center of the disk, but the difference becomes indistin- 
guishable in the spectrum of Procyon. 

3.2. Line asymmetries 

We have tried to avoid subjective criteria in selecting 
lines for measuring bisectors. Knowing that average con- 
vective line asymmetries for main-sequence F-K stars have 
been neither measured nor predicted to be larger than a 
few hundreds of meters per second, we have kept only those 
lines whose bisectors do not exceed 1 km s _1 in amplitude. 
The very large number of available lines in the spectrum 
whose rest wavelength is reliably known, makes it possi- 
ble to use such a strong rejection condition. In addition, 
only lines whose absorption was deeper than 10 % of the 
continuum level were included in the analysis. These cri- 
teria were satisfied by 331 lines in the solar atlas and by 
151 in Procyon's spectra. They were ordered in four dif- 
ferent groups depending on their continuum normalized 
flux at the line center h: h < 0.35, 0.35 < h < 0.55, 
0.55 < h < 0.75, and h > 0.75. Their bisectors were av- 
eraged retaining the relative velocity differences between 
the lowest end of the individual bisectors. In the process 
of averaging we allow for a cleaning, in the sense that at 
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each absorption depth at which the bisectors were mea- 
sured, we reject those deviating by more than 2a from the 
mean. 

Fig. H shows the averaged bisectors for Fe I lines. The 
solar bisectors are similar to the averaged line bisectors 
measured in the spectrum of the center of the disk by 
Dravins, Lindegren & Nordlund (1981). The typical 'C 
shape is clearly visible. A comparison of the velocity spans 
for the Sun and Procyon is in agreement with previous re- 
sults: the span of Procyon's lines exceeds the solar values 
by more than a factor of two. A similar result holds for the 
line shifts. As explained above, the velocity differences be- 
tween the individual bisectors were retained for averaging, 
but we have not subtracted the gravitational redshift in 
the solar case and the zero velocity for Procyon's bisector 
was arbitrarily set. 

The bisectors measured for Fe II lines are displayed in 
Fig. ^. As they are only a few, we did not average them 
out. Three solar bisectors (not shown) were rejected by the 
2a criterion, and there are still two more obvious outliers 
kept. Fe II bisectors corroborate the qualitative conclu- 
sions from Fe I lines. 

4. MODELING 

4.1. Fundamental stellar parameters 

Procyon A is in a 40-year period visual binary system, 
sharing the systemic velocity with a white dwarf. Gi- 
rard et al. (2000) have updated its astrometric orbit from 
photographic plates obtained at different observatories be- 
tween 1912 and 1995, and made a direct measurement 
of the angular separation of the pair using the infrared 
cold coronagraph (CoCo) at the NASA Infrared Telescope 
Facility and the WFPC2 onboard HST. At a distance of 
only 3.5 pc from the Sun, the parallax measured by Hip- 
parcos is very precise: p = (285.93 ± 0.88) x 10~ 3 arc- 
sec. In addition, Mozurkewich et al. (1991) measured 
the stellar angular diameter using optical interferometry: 



6 = (5.51 ± 0.05) 
be constrained to 



10 arcsec. Therefore, its mass can 



M 



a 



F 



,p2 



1.42±0.06M 



0! 



(2) 



where a, a a, and P are the semi-major axis of the vi- 
sual orbit or Procyon A, the semi-major axis of its or- 
bit relative to the barycenter of the system, and the or- 
bital period, respectively, and for which we have adopted 
a = 4.271 ± 0.032 arcsec, a A = 1-232 ± 0.08 arcsec, and 
P = 40.82 ± 0.06 yr. (Girard et al. 2000), in combination 
with Hipparcos^s parallax. The radius is 



R 



2ptan(^f) 



2.071 ± 0.020^ 



(3) 



where we have adopted 9q/2 = 959.64±0.02 arcsec (Choi- 
let & Sinceac 1999), leading to a gravity log g = 3.96±0.02 
(c.g.s. units). The derived gravity is in good agreement 
with the estimate of Allende Prieto & Lambert (1999), 
log<? = 4.04 ± 0.14 dex, obtained from matching the ob- 
served My and (B — V) to the evolutionary models of 
Bertelli et al. (1994) 5 . Our mass estimate is slightly lower 



from that given in Girard et al. (2000), as we prefer to 
adopt the Hipparcos parallax and the angular separation 
determined from HST-WFPC2 observations, discarding 
the ground-based measurement. Adopting the parallax 
derived by Girard et al., p = (283.2 ± 1.5) x 10~ 3 arc- 
sec, we find only a minor correction to the derived mass: 
M= 1.47±0.06M Q . 

The star has a Johnson visual magnitude of V = 0.363 ± 
0.003 mag (and therefore an absolute magnitude My = 
2.644 ± 0.007) and a color index B — V = 0.421 ± 0.003 
mag, as derived from 13 measurements extracted from 
SIMBAD. Fuhrmann et al. (1997) compiled several es- 
timates of Procyon's bolometric flux from the literature. 
Their weighted average is F BO l = (18.20 ± 0.43) x 10~ 6 



erg cm 
perature 



and, therefore, the star has an effective tern- 



Teg = 7.400 x 10 d 



F 



BOL 



1/4 



= 6512 ±49 K. 



(4) 



Finally, a review of the high-resolution spectroscopic anal- 
yses previously performed, indicates that the star has ei- 
ther solar, or slightly lower than solar, iron abundance. 

Linear interpolation in the solar-metallicity isochrones 
published by Bertelli et al. (1994) to match the pair 
(My,B — V) constrains very tightly the age of Procyon 
A to the range 1.66 — 1.73 Gyr, and predicts a gravity 
in agreement with the astrometric value. Alternatively, 
matching the observed radius and mass would shift the age 
estimate to 2.0 Gyr. This last value, less affected by un- 
certainties in the conversion between the theoretical and 
observational quantities, is to be preferred. These ages 
are significantly shorter than the nuclear time-scale of 3.9 
Gyr adopted by Provencal et al. (1997), implying a mass 
for the progenitor of the white dwarf (Procyon B) that is 
larger than previously estimated. 

4.2. Model atmospheres 

Realistic ab-initio, compressible, radiative- 
hydrodynamical simulations of surface convection in Pro- 
cyon have been performed with the code formerly used to 
model solar (e.g. Stein & Nordlund 1998; Asplund et al. 
2000a, b,c; Asplund 2000) and stellar granulation (Asplund 
et al. 1999; Asplund & Garcia Perez 2001). The result- 
ing structure is a 3D, time-dependent model atmosphere 
with a self-consistent description of the convective flows. 
The hydrodynamical equations of mass, momentum and 
energy conservation: 



<91np 
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coupled to the equation of radiative transfer: 
(n ■ V) I\ = r/ A - 



(5) 



(6) 



(7) 



(8) 



5 We note, however, that the estimate in a previous paper by Allende Prieto et al. (1999b) for this and other stars in their sample with 
T c ff > 6000 K is underestimated, as a result of the improper use of a different set of isochrones restricted to ages larger than 5 Gyr. 
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have been solved on a non-staggered Eulerian mesh with 
100 x 100 x 82 grid-points. In the above equations, p de- 
notes the density, v the velocity, g the gravitational ac- 
celeration, P the pressure, e the internal energy, a the 
viscous stress tensor, <5 v i sc the viscous dissipation, Q ra d = 
f. / n — S\)d£ldA the radiative heating/cooling rate, 

I\ the monochromatic intensity, with k\ and rj = B\(T)k\ 
representing the absorption and emission coefficients, re- 
spectively, where B\{T) is the Planck function, and n the 
unit vector for each of the eight considered directions. The 
code was stabilized using a hyper-viscosity diffusion algo- 
rithm (Stein & Nordlund 1998) with the parameters de- 
termined from standard hydrodynamical test cases, like 
the shock tube. Periodic horizontal boundary conditions 
and open transmitting top and bottom boundaries were 
used in an identical fashion to the simulations described 
by Stein & Nordlund. 

In order to make the atmospheric structure as realistic 
as possible for direct confrontation with observations, we 



used a detailed equation-of-state (Mihalas et al. 1988), 
including the effects of ionization, excitation and disso- 
ciation, as well as continuum (Uppsala opacity package; 
see Gustafsson et al. 1975; Asplund et al. 1997) and 
line (Kurucz 1993) opacities. The 3D radiative trans- 
fer has been solved at each time-step of the simulation 
by Feautrier's method, under the approximations of local 
thermodynamic equilibrium (LTE) and the opacity bin- 
ning technique (Nordlund 1982). The accuracy of the 
opacity binning procedure has been verified at regular in- 
tervals by solving the full monochromatic radiative trans- 
fer (2748 wavelength points) in the 1.5D approximation, 
i.e., treating each vertical column separately and thus ig- 
noring horizontal radiative transfer effects. 

The physical dimension of the numerical box is 21 x 21 x 
13 Mm, of which about 3 Mm correspond to the photo- 
sphere. A snapshot from a previous lower-resolution sim- 
ulation sequence of Procyon with a less extended depth- 
scale (Trampedach 1997) was used as initial condition. 
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The full convection simulation covers more than 3 hours 
of stellar time, but only the final hour has been used for 
the calculations of spectral line formation presented here. 
The resulting effective temperature of the 1 hour sequence 
was 6514 ±27K, in close agreement with the interferomet- 
ric estimate. We adopted a surface gravity of logy = 3.96 
(cgs) and solar chemical composition (Grevesse & Sauval 
1998). In particular, the Fe abundance used for the 
equation-of-state and continuum opacity calculations was 
log £fo = 7.50 6 . The line opacities were extrapolated to the 
adopted Fe abundance using the standard Kurucz ODFs 
with logepc = 7.67 and non-standard ODFs computed 
with logepc = 7.51 and with no He (Kurucz 1997, pri- 
vate communication, cf. Trampedach 1997 for details of 
the procedure). 

The resulting convective structures for Procyon are 
qualitatively similar to the lower resolution runs presented 
in Nordhmd & Dravins (1990), which, however, were com- 
puted with the anelastic approximation and less complete 
equation of state and opacities. With the new compress- 
ible code, larger convective velocities are encountered, and 
they often exceed the speed of sound, causing prominent 
shocks in the photosphere. The phenomenon of "naked 
granulation" discussed by Nordlund & Dravins (1990) is 
also present in the improved simulations as seen in Fig. 
[7]: the layers with largest temperature contrast due to 
the convective motions are located around the continuum 
forming region (rsooo — 1 — 3) and not slightly beneath it, 
as is the case with cooler stars like the Sun (t 50 oo — 10). 
In addition, the temperature contrast itself is significantly 
larger in Procyon than in the Sun. All these features 
manifest themselves in more pronounced line asymmetries. 
A detailed description of the convective structures in the 
granulation simulations of Procyon and other solar-type 
stars will be presented in a subsequent article, while here 
we concentrate on the effects upon spectral line formation 
and the confrontation with observations. 

To enable a strictly differential comparison between the 
3D and ID predictions of the resulting line formation, ID, 
hydrostatic, plane-parallel model atmospheres have been 
generated with the marcs code (Asplund et al. 1997). 
These homogeneous models are fully line-blanketed and 
have identical defining parameters as the 3D simulations, 
including all chemical abundances. 

4.3. Spectral line calculations 

Using the granulation snapshots of the Procyon simula- 
tion as 3D model atmospheres, the spectral line formation 
was computed for a sample of Fe I and Fe II lines assum- 
ing LTE (S v = B v ). The original simulation data was 
interpolated to a denser vertical grid covering only the 
outermost 6 Mm but with same number of depth points as 
before. From the stored temperatures and densities of the 
snapshots, monochromatic continuous opacities were com- 
puted using the Uppsala opacity package, which was also 
employed for constructing the ID model atmospheres we 
compare with. The assumptions of Boltzmann and Saha 
distributions and chemical equilibrium have been made in 
the calculation of opacities and number densities of the 
various species. 

The radiative transfer was solved for five center-to-limb 

6 On the customary logarithmic abundance scale defined to have log 



positions and seven equally spaced azimuthal angles, for a 
total of 29 different directions. The Gaussian distributed 
/i-angles were chosen to facilitate the disk integration tak- 
ing into account the rotational velocity of the star (Dravins 
& Nordlund 1990). 

The flux profiles at different velocities were computed 
as: 



F( v ) = / / H v ~ v rot sin i sin cosy;,//, y>)/id/idy>, 

(9) 

with fi = cos 6. The integrals were evaluated using a four- 
point Gaussian quadrature for \i and seven equidistant tp 
angles with equal weight. Finally, the rotationally broad- 
ened profiles for the 31 individual snapshots were added 
together before normalization to the local continuum. It 
was verified through test calculations that using twice as 
many fi- and y-angles results in insignificant differences 
in the theoretical line asymmetries. The procedure ac- 
counts for the Doppler shifts introduced by the convective 
motions in the stellar atmosphere which, together with 
the temperature inhomogeneities, produce the character- 
istic line bisector shapes. Various tests showed that the 
temporal coverage of the convection simulation produces 
statistically meaningful line profiles and asymmetries. Fol- 
lowing the calculations of disk-integration and rotational 
broadening, the resulting profiles were convolved with a 
Gaussian with a FWHM of 1.5kms _1 , to account for the 
finite spectral resolution of the spectrograph. As shown 
for the Sun, for which the rotational velocity is accurately 
known, there is no need to introduce any additional line 
broadening in the form of micro- and macroturbulence in 
analyses based on high-resolution 3D model atmospheres 
(Asplund et al. 2000a,b,c; Asplund 2000). 

The profiles were computed for 141 equidistant veloc- 
ities, covering —28 to +28kms~ 1 around the laboratory 
wavelengths. Each line was calculated for three different 
Fe abundances, 7.20, 7.50 and 7.80, from which the final 
profile was interpolated. Test calculations verified that the 
abundance interval was sufficiently small not to introduce 
any systematic errors in the abundances (< 0.01 dex) and 
line asymmetries (< 0.02 km s" 1 ) derived. 

The ID spectral line calculations for the comparison 
with the 3D results have been computed with the same 
spectral synthesis code as the 3D profiles. Without the 
convective Doppler shifts, additional adhoc broadening in 
the form of the micro- and macroturbulence must be in- 
voked in ID in order to obtain correct line widths, with 
the former affecting the line strengths and the latter only 
the line shapes. In both cases, Gaussian distributions are 
assumed. 

4.4. Line data 

Iron is the best represented element in the spectrum of 
late-type stars. Neutral iron has been the subject of a 
number of fine laboratory works to derive radiative tran- 
sition probabilities at Oxford (e.g. Blackwell et al. 1986). 
The number of lines measured with high accuracy and in 
a homogeneous manner has been enlarged by the work of 
O'Brian et al. (1991). Shown to be in good agreement 

<e h = 12.00 
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in Procyon. For each iso-tau surface, the temperature contrast (T rma /(T)) is given. 
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with the Oxford scale (see Lambert et al. 1996), we have 
adopted this last reference. Unfortunately, the situation 
for ionized iron is worse. To select transition probabili- 
ties for this species, we have restricted the search to lab- 
oratory measurements. Theoretical calculations still show 
many systematic effects as well as a large scatter. Astro- 
physical determinations always imply a certain prejudice 
that we wish to avoid: a selection of a model atmosphere, 
abundances, collisional damping parameters, etc. As we 
felt that an updated critical compilation of laboratory Fe 
II /-values was lacking, we have looked for data for all 
lines with wavelengths longer than 300 nm, rather than 
restricting the search to the spectral range of Procyon's 
observations. The method adopted to average the Fe II 
transition probabilities and the data are described in the 
Appendix. By combining laboratory, astrophysical, and 
theoretical determinations, it is possible to amass a larger 
line list, but at the expense of higher uncertainties (see, 
e.g., Giridhar & Arellano Ferro 1995). 

Whenever possible, new quantum-mechanical calcula- 
tions (Anstee & O'Mara 1991; Barklem & O'Mara 1997; 
Barklem et al. 1998, 2000) have been employed for the 
pressure broadening of the spectral lines by collisions with 
neutral hydrogen. Otherwise (in particular for the Fe II 
lines as the quantum-mechanical treatment has not yet 
been generalized to ionized species), the classical formula 
proposed by Unsold (1955), with an additional enhance- 
ment factor E = 2.0, has been adopted. Radiative damp- 
ing was included with values obtained from VALD (Kupka 
et al. 1999). Stark broadening was not considered. 

We selected lines of neutral and singly-ionized iron from 
the observed spectrum that were included in the list of 
O 'Brian et al. (Fe I), or in our list in the Appendix (Fe II). 
Many lines were eliminated from a preliminary list, since 
they were apparently blended with other features. Visual 
inspection decided which lines and parts of the profiles 
were selected at this point. The lines are listed in Table 1. 

In addition, we extracted the profiles of several Fe I and 
Fe II lines from the solar atlas of Kurucz et al. (1984), 
in order to estimate the solar iron abundance in a similar 
fashion. The lines were selected from the list in Asplund et 
al. (2000c). They show clean profiles, and the sources for 
the atomic data are identical to those used for Procyon. 
These lines and their atomic data appear in Table 2. In 
some cases, the transition probabilities differ from those 
given in Asplund et al. (2000c); the difference is however 
marginal (~ 0.02 dex on average for both Fe I and Fe II). 
We noticed some solar Fe lines partially affected by blends, 
and the afflicted parts of the profiles were ignored in the 
analysis. 

5. COMPARISON BETWEEN THEORETICAL AND 
OBSERVED LINE PROFILES 

5.1. Line profiles and stellar rotation 

In the case of 3D analyses, u ro tsim can be determined by 
minimizing the residuals between theoretical and observed 
profiles, since no other unknown broadening like macro- 
and microturbulence enters the calculations. This is per- 
formed through a x 2 -analysis with the elemental abun- 
dance, ^rotsini, FWHM, and t> ra d as free parameters. By 

7 For the FTS solar observations the instrumental profile is better 
calculated profiles. 



FWHM we refer to that of an assumed-Gaussian profile 
that in 3D accounts for the instrumental profile, while w ra d 
is the required overall velocity shift of the whole observed 
profile. In the similar ID analysis described below the 
FWHM represents the combined effect of instrumental res- 
olution and macroturbulcncc, while the microturbulence is 
kept fixed in the x 2 -analysis. Figs. || and [)] compare the 
observations and the best-fit synthetic spectra calculated 
with ID and 3D model atmospheres, respectively, for sev- 
eral iron lines. A common deficiency of the ID predictions 
is apparent from Fig. 8, where the residuals show an os- 
cillation around the line center. Fig. [l(] shows the best-fit 
parameters obtained from individual lines. 

For Procyon, we estimate D ro tsmi = 3.16 ± 0.50 km s , 
where the quoted uncertainty is the standard deviation 
for the clean samples of Fe I and Fe II lines minus a few 
obviously discrepant cases (see Fig. |l0|). There is, except 
for a couple of aberrant lines, no need for additional Gaus- 
sian broadening above the known instrumental broadening 
(1.5kms _1 ); macroturbulence, which must be introduced 
with the ID model, is not required. The residual fluxes 
show no systematic behavior across the profile, as the cor- 
responding ID lines do. In fact, the 3D modeling makes it 
possible to detect blended lines which would have escaped 
detection in ID analyses. 

It should be noted, however, that the above estimate of 
^rotsim may be slightly over-estimated due to the finite 
numerical resolution of the convection simulation. With 
higher resolution more of the high- velocity tails of the ve- 
locity distributions are sampled, which causes additional 
broadening of the lines (Asplund et al. 2000a). For the 
Sun, the predicted line profiles appear to have converged at 
the currently highest affordable resolution (200 x 200 x 82). 
Since the present Procyon simulations are based on a grid 
size of 100 x 100 x 82 due to computing time considera- 
tions, the widths of the theoretical line profiles may have 
been slightly under-estimated. Such a conclusion is sup- 
ported by the presence of a persistent trend in derived Fe 
abundances with line strength (Sect. |^). Assuming a sim- 
ilar behavior with numerical resolution for Procyon as for 
the Sun, w rot sim may be « 0.5 km s -1 less than estimated 
above. However, this additional broadening has a very lim- 
ited impact on the predicted line shifts and asymmetries 
(Asplund et al. 2000a), in particular for weak lines, which 
are still reliable indicators of the chemical abundances. 

A detailed Fourier analysis by Gray (1981b), using 
classical ID model atmospheres and McDonald spectra 
with a resolution of 2.5 km s -1 , provided w ro tsim = 2.8 ± 
O.Skms" 1 , and a radial-tangential macroturbulence of 
7.0 ± 0.1 kms" 1 . Studies by Fekel (1997) and Benz & 
Mayor (1984) gave estimates for u rot sini of 4.9 ± 1.0 and 
4.5 ± 1.1 kms -1 , respectively. Our own x 2 -analysis us- 
ing ID model atmospheres reveals similar values but also 
that the derived i> ro tsini depends on the choice of line, in 
particular for weak lines. At least partly, this will de- 
pend on the form of the adopted macroturbulence, for 
example, Gaussian as here, or radial-tangential as used 
by Gray (1981b). The scatter in w rot sim is clearly signif- 
icantly larger with ID than with 3D model atmospheres 
when using otherwise identical analyses. The reason why 
lines with W\ ~ 70 — 100 mA apparently indicate the cor- 

described by a sine function, but this has a negligible impact on the 
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Fig. 8. — Comparison between the synthetic profiles calculated with the ID model atmosphere of Procyon (solid line) against the 
observations (filled circles). The residuals (predicted - observed) are also shown, but multiplied by a factor of 5. 
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Fig. 9. — Comparison between the synthetic profiles calculated with the 3D model atmosphere of Procyon against the observations (filled 
circles). The residuals (predicted - observed) are also shown, but multiplied by a factor of 5. 
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rect ^rotsim, while weaker lines do not, remains a mystery 
to us; in fact the solar analysis shows the same effect (Fig. 
lit). 



5.2. Line asymmetries and shifts 

A major advantage with using 3D hydrodynamical 
model atmospheres for spectral synthesis besides the elim- 
ination of free parameters like mixing length parame- 
ters, micro- and macroturbulence, is that the convective 
velocity field and temperature inhomogeneities are self- 
consistently computed. As a consequence, the imprint of 
convection in the form of asymmetries and shifts of spec- 
tral lines can be predicted and compared with observa- 
tions. Furthermore, the predicted line asymmetries are 
on an absolute wavelength scale. Since different lines have 
different sensitivity to the photospheric structure and con- 
vective Doppler shifts, all lines show unique asymmetries. 

Figure [l| illustrates the excellent agreement between 
predicted and observed line asymmetries for both Fe I and 
Fe II lines and weak and strong lines. Since the exact 
radial velocity of Procyon is unknown, the observed spec- 
trum is not on an absolute wavelength scale. Therefore, 
the zero-point for all observed bisectors in Fig. [l2| has been 
shifted by the same amount to make them coincide with 
the theoretical bisectors on average (see further discussion 
below and in Sect. 8). It should be noted, however, that 
the observed relative shifts between different lines shown 
in Fig. O have been derived from the observations and 



are obviously well predicted by the 3D calculations. 

Although the general agreement between theory and ob- 
servations is quite satisfactory, there are several notable 
exceptions for which the discrepant bisectors may signal 
problems like blends or an erroneous laboratory wave- 
length. In the case of Fe I 4745.8 A, Fe I 5434.5 A, and 
Fe I 4576.3 A the reason is most likely minor blends in 
the wings, which is supported by the good correspondence 
between predicted and observed bisectors for larger line 
depths. The Fe II 4629.339 A line can be suspected to 
have an erroneous laboratory wavelength due to its un- 
characteristically large offset of about 9 mA to the red. 
It is noteworthy that this line is the only Fe II line in 
our sample which is not included in the compilation of S. 
Johansson (1998, private communication) and the labora- 
tory wavelength has instead been adopted from the VALD 
database. 

The bisector comparison is summarized in Figure [l3[ 
which shows the difference between the theoretical and ob- 
served bisectors for the clean samples of Fe I and Fe II lines 
together with a few additional unblended Fe II lines which 
lack an accurate determination of the transition probabil- 
ity and therefore are not included in the estimate of the 
Fe abundance presented below (see §^|). For most lines, 
the predicted bisector is simply offset from the observed 
bisector, i.e. the difference is essentially a vertical line, 
without trend with line depth that could signal an im- 
proper theoretical velocity field. The scatter around zero 
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difference is, at least partly, due to errors in the laboratory 
wavelengths, as well as the finite S/N and spectral reso- 
lution of the observations. Close to the continuum, the 
effects of weak blends become significant, which causes a 
larger scatter and in a few cases a deviating bisector. Typi- 
cally, the shapes of the bisectors can be predicted to within 
50ms" 1 , while the typical velocity spans are 500ms _1 or 
more, which is quite remarkable considering the sensitivity 
of line asymmetries upon the detailed atmospheric struc- 
ture. In fact, the agreement is even slightly better than 
for the Sun (Asplund et al. 2000b). 

A particularly interesting result arising from the line 
asymmetry comparison, is the predicted blueward hook 
close to the continuum in the bisectors of some lines. All 
predicted Fe lines with h < 0.5, and hence partly satu- 
rated, show this behavior, while it is not present in any 
weaker lines. The blue hook is induced by the influence of 
extended Lorentzian wings from the upflowing and hence 
blue-shifted granules (Dravins 1990). Due to the bias in 
line strength towards granules as a result of their steeper 
temperature gradients, lines tend to first saturate and de- 
velop damping wings there. The extended nature of these 
blue-shifted wings can thus influence the resulting bisector 
shape, in particular close to the continuum. Even though 
blends tend to increase the uncertainty in the measured 
bisectors close to the continuum, there also exists obser- 
vational evidence for this blue-ward hook. In particular 
Fe I 5415.2 A is a convincing case but also Fe I 4903.3 A, 
Fe I 5397.1 A, and Fe II 4620.5 A support this conclusion. 



From the Procyon atlas of Griffin & Griffin (1979), Dravins 
(1987) claimed to have identified a similar feature by aver- 
aging over many lines, however, the feature is n ot apparent 
in the average bisectors discussed in Sect. 3.2. 

As expected from the good general agreement in terms 
of line asymmetries, the predicted differential line shifts 
match closely to the observations. Fig. |lj shows the 
difference between predicted and observed line shifts for 
the Fe I and Fe II lines; the Fe I 4903.3101 A and Fe 
II 4629.3390 A lines stand out; probably with erroneous 
laboratory wavelengths. Since the radial velocity of Pro- 
cyon is unknown, the observed profiles have been shifted 
by — 1.56kms _1 in order to make the difference between 
predictions and observations disappear on average. The 
mean of the differences between predicted and observed 
line shifts amounts to O.OOkms -1 with a = O.lOkms" 1 , 
implying that the radial velocity can be equally well de- 
termined in 3D analyses from fitting the whole profiles 
(Fig. |hJ) or only the line center. The observed weak- 
ening of the convective blue-shift for stronger lines is the 
result of the the disappearance of the convective inhomo- 
geneities as the line-formation region is gradually moved 
outwards. Unfortunately, the present observations do not 
include sufficiently strong lines to determine if a plateau 
at zero convective line shifts exists, like that found for 
the Sun (Allende Prieto & Garcia Lopez 1998), and the 
3D convection simulation can not be invoked either to an- 
swer this question. As clear from Fig. the predicted 
line shifts become systematically biased for stronger lines, 
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which is not surprising given the likely departures from 
LTE in the cores of such lines. Besides, the 3D simula- 
tions lose realism for the upper atmosphere as a result of 
the influence of the upper boundary and missing magnetic 
fields. In those layers, the presence of propagating waves 
and the fact that the energy balance is determined by a 
few strong lines out of LTE complicates quite seriously the 
modeling. A similar deficiency was detected in the solar 
study (Asplund et al. 2000b). 

In Figure 15 we have mimicked Fe I lines of different 
strengths by altering the iron abundance in the calcula- 
tions. Although a quantitative comparison is not possible, 
because of the different ways in which the two plots were 
constructed, the general behavior of Fe I line asymmetries 
shown in Figure |5j is clearly reproduced by the theoretical 
calculations. We have calculated the line velocity shifts (v) 
for a large group of Fe I lines, and as Figure [l(] reveals, 
the intrinsic scatter of the convective shifts for weak lines 
about a straight fit is significantly smaller in the Sun than 



in Procyon. This r esult confirms qualitatively the obser- 
vational findings in £3.1. It should be noted, however, that 
since the same lines have not been used for the theoretical 
and observed distributions one should not be surprised of 
minor differences in the results. The predicted slopes of 
the velocity shifts as a function of equivalent width are 
5.9 x 10 -3 for the Sun and 6.0 x 10~ 3 for Procyon - vir- 
tually identical. Neither of them are very close to the 
observed slopes, but the agreement is clearly better for 
Procyon, whose observed value is 4.9 x 10~ 3 , than for the 
Sun (3.2 x 10~ 3 ; see Asplund et al. 2000a). 

The line asymmetry comparison clearly shows that the 
adopted 3D model atmosphere provides quite a realistic 
description of the granulation properties on Procyon. We 
stress again that ID models can predict neither the bisec- 
tors nor the shifts of lines. 



6. THE IRON ABUNDANCE OF PROCYON 
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The chemical abundances inferred from the comparison 
between our spectroscopic observations and the 3D hy- 
drodynamical model atmosphere of Procyon have a spe- 
cial interest, as this is the most common driver for stellar 
spectroscopy. It is particularly useful to refer our results 
to those from classical ID model atmospheres, on which 
most of the published stellar abundances are based. Here 
we determine the iron abu ndance of Procyon with the two 
types of models (see |4.2| ), making a strictly differential 
comparison between them. 

6.1. Classical ID models 

In addition to the free parameters in the x 2 -analysis 
(iron abundance, w ro tsinz, FWHM, and w ra d), one must 
also specify the microturbulcnce (£tur) in analyses based 
on ID model atmospheres. Variations in this parameter 
affect both the line strengths and shapes, in particular for 
partly saturated lines. Following standard practice, we 
have carried out several x 2_an alyses using different val- 
ues of £tur in order to select the value which removes 
any trend in derived abundances with line strength. For 
Procyon, our choice of lines, transition probabilities, and 
model atmospheres suggest a value of £ tur = 2.2 km s -1 . 



The results are displayed in the right-hand panels of Fig. 



17. The mean abundance derived from the neutral- iron 
lines (filled circles) and ionized-iron lines (open circles) is 
loge(Fe) = 7.30 ± 0.02 (a = 0.11) dex and 7.32 ± 0.03 
(a = 0.08) dex, respectively. 

A similar analysis of the same lines used by Asplund et 
al. (2000c) in the disk-integrated solar atlas (Kurucz et 
al. 1984), but adopting the same sources for the radiative 
transition probabilities used for Procyon's (§ [1.4| ), provided 
loge (Fe) = 7.47 ± 0.01(ct = 0.07) dex for the Fe I lines 
and loge (Fe) = 7.41 ± 0.03(er = 0.11) dex for the Fe II 
lines with a MARCS model atmosphere and £ = 1.0 km s _1 . 
These figures are represented against equivalent width and 
excitation potential in the right-hand panels of Fig. |l8|. 

6.2. Three-dimensional models 

Without an adjustable microturbulence parameter, the 
process of determining the abundances from the line pro- 
files is straightforward from the x 2 -analysis. To take ad- 
vantage of all the information in the line profiles, and in 
parallel to our procedure for the ID models, we fit v rot sin i, 
FWHM, and u rad , for each individual line. Now, FWHM 
is also that of a Gaussian, but it only accounts for the in- 
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strumental prof ile, as macro-turbulence is deemed as un- 
necessary (§5.1). The left-hand panels of Figure [Tt] show 
the abundances from different lines plotted against the 
line equivalent width and the line excitation potential. 
The mean abundance for the neutral iron lines is log e 
(Fe) Q = 7.41 ± 0.02 {a = 0.14) dex, and for the singly- 
ionized iron lines loge (Fe) Q = 7.49 ± 0.04 (a = 0.10) 
dex. The retrieved abundances tend to be higher for the 
stronger lines, as also found for the solar case with a more 
limited range of equivalent widths (see Fig. |l8|) . This can 
be likely ascribed to the limited resolution of the simu- 
lations, which miss the highest convective velocities and 
therefore Doppler shifts. To put this into perspective we 
note that this trend is similar to an underestimated mi- 
croturbulence of about 0.3 km s _1 for Procyon, which 
highlights that it is a relatively minor shortcoming of the 
3D simulations as most of the needed convective Doppler 
shifts are already accounted for. Restricting the average 
to lines with equivalent widths less or equal than 50 mA 
we obtain 7.35 ± 0.03 dex (a = 0.15 dex; 22 lines) and 
7.36 ± 0.01 dex (a = 0.01; 2 lines) for neutral and singly- 
ionized iron lines. 

The analysis of the selected solar flux lines previously 
used for the ID case, together with the 200 x 200 x 82 
3D simulation described by Asplund et al. (2000b), led 
to the following estimates of the solar iron abundance: 
loge (Fe) Q = 7.40 ± 0.01 (a = 0.07) and 7.43± 0.03 
(er = 0.11) dex, as derived from the neutral and singly- 
ionized iron lines, respectively. These figures change very 
mildly when we restrict the sample to lines with equiva- 
lent widths smaller than 50 mA, as all the lines considered 



have equivalent widths smaller than 100 mA. 

The small differences in derived abundances with those 
given in Asplund et al. (2000c), which are based on the 
same 3D model atmosphere as employed here, can be 
traced to slightly different adopted transition probabili- 
ties (<; 0.02 dex) and the use of flux profiles instead of 
disk-center intensity profiles (<; 0.02 dex). Our new esti- 
mates do not supersede the solar iron abundances given in 
Asplund et al. (2000c). They are only intended to enable 
a differential study of Procyon. 

6.3. Discussion 

Estimates of the accuracy of the log gf for the transi- 
tions used are readily available, either from O 'Brian et al. 
(1991), or from the tables in the Appendix. We can use 



them to perform an analysis similar to that in §3.1 for 



the line velocity shifts. This reveals that the errors in the 
logg/s cannot account for more than about half of the 
scatter in the iron abundances. Even though 3D models 
need to be improved, there are probably other important 
error sources in our calculated line profiles. Departures 
from LTE are a likely candidate. We do not observe a 
larger scatter or a systematic deviation for low-excitation 
Fe I lines that one could expect from the 3D NLTE cal- 
culations of Trujillo Bueno & Shchukina (2001), but we 
notice a significantly larger scatter for the iron abundance 
determined from Fe I lines in Procyon than in the Sun, 
and departures from LTE are expected to grow for warmer 
stars. 

Because of the shown imperfections in the 3D mod- 
els, the preferred iron abundance for Procyon is that de- 
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Fig. 17. — Iron abundance derived for Procyon from Fe I (filled circles) and Fe II lines (open circles), as a function of the equivalent width 
(upper panel) and the line excitation potential (lower panel). 



rived from the 3D analysis of weak lines: loge (Fe)© = 
7.36 ±0.03 (ct = 0.15). 

The main conclusions of our study of the iron abundance 
are: 

1 . The abundances derived from Fe I and Fe II lines are 
consistent both for ID and 3D analyses, and with 
a higher coherence for the 3D model. This implies 
that departures from LTE, which should be minimal 
for Fe II, are likely small. Nevertheless, they could 
be responsible for a significant fraction of the scat- 
ter between the abundances retrieved from different 
lines, providing a plausible explanation for the larger 
scatter observed for Procyon than for the Sun. 

2. The differences between the iron abundances derived 
from ID and 3D analyses (setting aside lines stronger 
than 50 mA in the case of Procyon) are small (<J 0.05 
dex), and the same conclusion applies to the solar 
case. Other stars and elements may show larger dif- 
ferences. 

3. Procyon is marginally deficient in iron compared to 
the Sun by about 0.05 dex. 

7. CENTER-TO-LIMB VARIATION 

One of the applications of model atmospheres is to de- 
rive limb-darkening laws. These are required, for example, 



to correct interferometric measurements of stellar angular 
diameters, and affect directly the otherwise fully-empirical 
calibrations of effective temperature against color indices 
(see, e.g. Mozurkewich et al 1991). 

Here we compare the center-to-limb variation predicted 
by homogeneous models and the 3D simulations. Taking 
the emerging intensity for different angles at all different 
spatial locations we produce a spatially and time-averaged 
intensity for different ray inclinations. Figure [l^ shows the 
predicted limb-darkening for the ID and 3D models. Two 
wavelengths were selected, 4500 and 10000 A. Strong dif- 
ferences are obvious in the plot, the limb-darkening being 
markedly non-linear for the inhomogeneous model. We 
have fitted the data to third-order polynomials by regu- 
lar least squares (see Equation 6 in Hanbury Brown et 
al. 1974) and the results correspond to the solid (3D) 
and dashed (ID) lines in Figure [n]. Lower order polyno- 
mials, were inappropriate for the 3D model atmosphere. 
Hanbury Brown et al. (1974) calculated the interferomet- 
ric correlation factor associated with this particular limb- 
darkening law. At 4500 A, their determined angular diam- 
eter for a uniformly emitting disk should be corrected by 
factors of 1.081 and 1.064 for the ID and 3D cases, respec- 
tively. The ID correction obtained from an older version 
of MARCS (Gustafsson et al. 1975) used by Mozurkewich 
et al. (1991) was 1.07 at 4500 A. The radius and the ef- 
fective temperature of Procyon we derived in §4.1 should 
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therefore be slightly corrected from 2.071 to 2.059 Rq and 
from 6512 to 6530 K, respectively. The correction between 
3D and ID, which amounts to roughly 1.6 %, implies a cor- 
rection to the effective temperature of roughly 50 K for a 
star like Procyon A, showing the importance of detailed 
model atmospheres for establishing a truly empirical T c g 
scale from measurements of angular diameters (see, e.g., 
Di Benedetto 1998). 

8. ABSOLUTE RADIAL VELOCITY 

When discussing measurements of stellar radial veloc- 
ities, it is commonly assumed that an observed shift of 
the spectral features from their wavelengths at rest can be 
directly translated to a velocity projected along the line 
of sight. Lindegren, Dravins & Madsen (1999) have re- 
cently discussed the dangers in doing so, especially when 
it comes to precise measurements, pointing out effects that 
may need to be considered. 

The gravitational shift is proportional to the gravita- 
tional field, V g = GM/(Rc), and for late-type dwarfs 
amounts to about 0.5 km s _1 . In the case of Procyon 
A, we have a special advantage, as its mass, radius, and 
proper motion are well determined. The gravitational red- 
shift, including the blueshift induced by Earth is V g = 
0.436 ± 0.019 km s , where we have used the corrected 
stellar radius R = 2.06 ± 0.02i? Q (see §0). The transverse 
Doppler effect scales with the square of the modulus of the 



velocity, Vto — | | 2 / (2c) , and becomes significant, reach- 
ing ~ 0.25 km s _1 for the fastest halo stars orbiting at ~ 
400 km s . The proper motion of Procyon as measured 
by Hipparcos corresponds to a velocity perpendicular to 
the line of sight of Vt = 20.87 ± 0.07 km s , which trans- 
lates to a negligible transverse Doppler shift of < 1 m s . 

The list of contributors to the observed line shifts that 
have little to do with the star's space motion is much 
longer: pulsations, stellar winds, surface convection, pres- 
sure gradients, etc. In the case of quiet late-type dwarfs, 
stellar winds are not expected to disturb the photospheric 
line wavelengths, and the integrated light pulsations are 
known to have periods of minutes and amplitudes of a few 
meters per second or less, as detected in the Sun (see, e.g., 
Fossat & Ricort 1975, Claverie et al. 1979) or Procyon 
itself (Martic et al. 1999). Pressure gradients have been 
shown to produce line asymmetries and shifts of only a 
few meters per second for solar-like stars (Allende Prieto, 
Garcia Lopez & Trujillo Bueno 1997). In conclusion, the 
main disturbing element is the velocity pattern of convec- 
tive origin. The detailed three-dimensional hydrodynam- 
ical simulations here presented provide an accurate de- 
scription of the convective velocity-temperature fields in 
the photosphere, and therefore constitute a means to zero 
the photospheric velocities. Confronting the measured line 
shifts with the predictions of the 3D model (see Section 
5.2), and restricting the analysis to lines weaker than 100 
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Fig. 18. — Iron abundance derived for the Sun from Fe I (filled circles) and Fe II lines (open circles), as a function of the equivalent width 
(upper panel) and the line excitation potential (lower panel). 
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mA, we determine a velocity shift of the observed spec- 
trum of Procyon A by —1.56 ±0.13 km s _1 . That amount 
has to be divided between the remaining effects expected 
to be important in the case of Procyon A: gravitational 
rcdshift, orbital velocity, and systemic velocity. 

Combining the astromctric orbital elements of Girard 
et al. (1996) with the radial velocity amplitude measured 
by Irwin et al. (1992): 1.401 ± 0.093 km s" 1 (from abso- 
lute measurements) or 1.482 ±0.110 km s _1 (from relative 
determinations), we find that the orbital velocity of Pro- 
cyon A changed by less than 0.0004 km s _1 during the 
two days the observations were carried out, and its value 
was 1.115 (from absolute measurements) or 1.180 km s _1 
(from relative determinations). Adopting 1.15 ± 0.10 km 
s _1 , we obtain for the heliocentric radial velocity of the 
system: V r = -1.56(±0.13) - 0.44(±0.02) - 1.15(±0.10) = 
-3.15 ±0.17 km s" 1 . 

9. SUMMARY AND CONCLUSIONS 

We have observed the spectrum of Procyon A (F5IV) 
from 4559 to 5780 A with a S/N of ~ 10 3 and a resolving 
power of 2 x 10 5 . We have measured the line bisectors and 
relative line shifts of large number of Fe I and Fe II lines. 
Compared to the Sun, the convective velocity shifts of the 
lines arc significantly larger and show a larger scatter from 
a linear relationship against equivalent width. The differ- 
ence between the weakest and the strongest lines measured 
can reach up to 1 km s _1 , which is almost twice as large as 
the range observed in the Sun. Among other interesting 
consequences, this large differential velocities are likely to 
limit the accuracy in the determination of absolute radial 
velocities in F-type stars. 

We have constructed a new three-dimensional, time- 
dependent hydrodynamical model atmosphere, and tested 
it against the observed spectrum. The new model repro- 
duces in detail most of the observed features, but we iden- 
tify some room for improvement. At all levels, the com- 
parison with the observed spectral lines shows a significant 



improvement compared to classical homogeneous models. 
The predicted line shifts and asymmetries are reasonably 
close to the observed ones. Our study of the iron abun- 
dance in the photosphere of Procyon reveals that: 

• The abundances derived from Fe I and Fe II lines 
are consistent both for ID and 3D analyses, and to 
a higher degree for the latter. This implies that de- 
partures from LTE, which should be minimal for Fe 
II, are relatively small. 

• the differences between the mean iron abundances 
derived from ID and 3D analysis are small (<; 0.05 
dex). 

• we find the iron abundance of Procyon to be 
loge(Fc) = 7.36 ±0.03 {a = 0.15) dex, concluding, 
from a differential study of the solar spectrum, that 
the star is marginally deficient in iron compared to 
the Sun by ~ 0.05 dex. 

We notice a scatter in the abundances from different 
lines that clearly exceeds the expectations from errors in 
the transition probabilities. Known imperfections in the 
3D model should account for part of the excess errors. De- 
partures from LTE are a potentially important source of 
error. 3D models with improved resolution and detailed 
NLTE calculations are needed. In addition, other deficien- 
cies in the input physics could contribute, and could be 
identified through meticulous tests against observations. 

The detailed line shapes predicted by the model provide 
a reliable upper limit for the small projected rotational 
velocity of the star, v Iot sinz < 3.16 km s _1 , with the cor- 
rect value probably close to 2.7 km s _1 . The 3D model 
atmosphere shows a limb-darkening law that is noticeably 
different from the predictions of ID models. This finding 
has limited impact on the interferometric determination 
of the angular diameter of Procyon (~ 2 %) and its ef- 
fective temperature (~ 50 K), but is obviously important 



in establishing an absolute scale of T c g based on direct 
measurements. The ability of the hydrodynamical model 
atmosphere to predict the convective shifts of the spec- 
tral lines, combined with the accurate parallax and proper 
motions determined by Hipparcos, the published orbital el- 
ements, and a careful assessment of the gravitational red- 
shift makes it possible to determine the space motion of 
the Procyon binary system within 0.17 km s _1 . 
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APPENDIX 

A Compilation of Laboratory Absolute Transition Probabilities for Fe II lines with A > 300 ran. 

The experimental technique most used combines measured branching ratios and lifetimes of the upper levels to derive 
transition probabilities. Measurements of branching ratios of high quality have been published by Kroll & Kock (1987), 
Whaling (as included in the compilation by Fuhr, Martin & Wiese 1988 - also private communication), Heise & Kock 
(1990), and Pauls, Grevesse, & Huber (1990). Level lifetimes have been derived with accuracies in the range 0.02-0.14 
ns by Biemont et al. (1991), Guo et al. (1992), and Schnabel, Kock & Holweger (1999) using laser-induced fluorescence. 
Available lifetimes have been compiled and combined in Table |^. When more than one measurement was available, we 
averaged them out. The variation among the lifetime of the levels within a term is remarkably small (except maybe 
for the 4 F term), as expected if LS-coupling holds. Lower accuracy measurements that include 4 Di/ 2 and 4 F 3 / 2 levels 
(providing for them 2.9 ± 0.2 and 3.7 ± 0.2, respectively), support the use of the term average values for those cases 
(2.95 ± 0.03 and 3.70 ± 0.15). The branching fractions measured by Kroll & Kock, Whaling, Heise, & Kock, and Pauls 
et al. were converted to transition probabilities using different sources for the lifetimes of the involved upper levels. The 
third through sixth columns of Table |J list their published log(g/)s, the quoted uncertainties, and the lifetimes used for 
the different transitions. We have re-scaled the individual log(g/)s using the adopted lifetimes of Table ^, to produce a 
weighted averaged when more than a single measure was available. When inconsistent log(<?/)s were derived from Kroll 
& Kock and Heise & Kock, the latter were preferred, following Schnabel et al. The finally adopted log(<?/)s and they 
formal uncertainties are displayed in the last column of Table 0. 
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Table 1 

Spectral Lines Selected for the Iron Abundance Analysis of Procyon 



Species 


Wavelength 


ilir 


log gj 


log e oJJ 


loge 1JJ 


Wa(3L>) 




(A) 


(eV) 


(dex) 


(dex) 


(dex) 


(mA) 


Fe I 


4602 00 


1.608 


-3.13 


7.30 


7.31 


47 


Fe I 


4602 94 


1.485 


-2.21 


7.51 


7.27 


105 


Fe I 


4635 85 


2.845 


-2.36 


7.32 


7.34 


31 


Fe I 


4683 56 


2.831 


-2.32 


7.24 


7.25 


30 


Fe I 


4690.14 


3.686 


-1.64 


7.49 


7.44 


39 


Fe I 


4710.28 


3.018 


-1.61 


7.52 


7.41 


72 


Fe I 


4735.84 


4.076 


-1.33 


7.65 


7.57 


46 


Fe I 


4741.53 


2.831 


-1.76 


7.16 


7.11 


54 


Fe I 


4745.80 


3.654 


-1.27 


7.47 


7.39 


58 


Fe I 


4779.44 


3.415 


-2.02 


7.22 


7.23 


21 


Fe I 


4786.81 


3.017 


-1.61 


7.50 


7.39 


70 


Fe I 


4788.76 


3.237 


-1.76 


7.39 


7.35 


48 


Fe I 


4789 65 


3.546 


-0.96 


7.49 


7.31 


77 


Fe I 


4798.27 


4.186 


-1.17 


7.16 


7.12 


25 


Fe I 


4802 88 


3.642 


-1.51 


7.40 


7.36 


43 


Fe I 


4832 73 


3.640 


-1.73 


7.61 


7.57 


42 


Fe I 


4844.01 


3.546 


-2.05 


7.46 


7.45 


25 


Fe I 


4872.14 


2.882 


-0.57 


7.45 


7.22 


133 


Fe I 


4903 31 


2.882 


-0.93 


7.64 


7.40 


120 


Fe I 


4918 99 


2.865 


-0.34 


7.52 


7.33 


160 


Fe I 


4985 25 


3.928 


-0.56 


7.49 


7.27 


85 


Fe I 


4994.13 


0.915 


-2.97 


7.26 


7.15 


87 


Fe I 


5006 1 2 


2.832 


-0.62 


7.61 


7.48 


143 


Fe I 


5014.94 


3.943 


-0.30 


7.48 


7.26 


98 


Fe I 


5022 79 


2.990 


-2.20 


7.46 


7.44 


41 


Fp T 


'J\J L ±iJ . 


9 979 


1 36 


7 59 


7 30 


1 1 8 

no 


Fe I 


5054.64 


3.640 


-1.92 


7.20 


7.21 


17 


Fe I 


5079.74 


0.990 


-3.24 


7.42 


7.31 


78 


Fe I 


5083.34 


0.958 


-2.84 


7.25 


7.12 


91 


Fe I 


5127.36 


0.915 


-3.25 


7.28 


7.23 


75 


Fe I 


5141.74 


2.424 


-2.24 


7.41 


7.34 


62 


Fe I 


5145.09 


2.198 


-2.88 


7.05 


7.08 


22 


Fe I 


5151.91 


1.011 


-3.32 


7.35 


7.30 


71 


Fe I 


5194.94 


1.557 


-2.02 


7.36 


7.14 


109 


Fe I 


5198.71 


2.223 


-2.09 


7.35 


7.24 


77 


Fe I 


5232.94 


2.940 


-0.06 


7.41 


7.25 


175 


Fe I 


5242.49 


3.634 


-0.97 


7.46 


7.31 


75 


Fe I 


5288.52 


3.694 


-1.51 


7.31 


7.29 


37 


Fe I 


5321.11 


4.434 


-1.09 


7.22 


7.20 


23 


Fe I 


5322.04 


2.279 


-2.80 


7.22 


7.26 


31 


Fe I 


5328.04 


0.915 


-1.47 


7.50 


7.29 


184 



Table 1 — Continued 



Species 


Wavelength 


EP 


log gf 


log e 3D 


loge ID 


W A (3D) 




(A) 


(eV) 


(dex) 


(dex) 


(dex) 


(niA) 


Fe I 


5364.87 


4.445 


0.23 


7.51 


7.31 


110 


Fe I 


5365.40 


3.573 


-1.02 


7.21 


7.11 


63 


Fe I 


5367.47 


4.415 


0.44 


7.43 


7.22 


121 


Fe I 


5383.37 


4.312 


0.64 


7.47 


7.24 


139 


re 1 


COOT 1 O 

5397.13 


0.915 


1 AO 

-1.9o 


( .60 


7 .29 


152 


Fe I 


5403.82 


4.076 


-1.03 


7.49 


7.38 


54 


Fe I 


5415.20 


4.386 


0.64 


7.40 


7.21 


138 


Fe I 


5434.52 


1.011 


-2.13 


7.59 


7.29 


139 


Fe I 


5446.92 


0.990 


-1.91 


7.64 


7.33 


157 


Fe I 


5464.28 


4.143 


-1.40 


7.21 


7.21 


21 


Fe I 


5466.99 


3.573 


-2.23 


7.55 


7.58 


21 


Fe I 


5483.10 


4.154 


-1.41 


7.37 


7.35 


27 


Fe I 


5497.52 


1.011 


-2.83 


7.55 


7.45 


105 


Fe I 


5501.47 


0.958 


-3.05 


7.48 


7.33 


94 


Fe I 


5586.76 


3.368 


-0.14 


7.53 


7.30 


144 


Fe I 


5618.63 


4.209 


-1.27 


7.38 


7.36 


31 


Fe I 


5701.54 


2.559 


-2.14 


7.34 


7.29 


59 


Fe I 


5741.85 


4.256 


-1.67 


7.41 


7.41 


16 


Fe II 


4576.33 


2.844 


-2.91 


7.52 


7.26 


98 


Fe II 


4620.51 


2.828 


-3.19 


7.42 


7.25 


80 


Fe II 


4629.34 


2.807 


-2.28 


7.54 


7.24 


133 


Fe II 


4923.92 


2.891 


-1.26 


7.50 


7.41 


234 


Fe II 


5234.62 


3.221 


-2.23 


7.62 


7.36 


127 


Fe II 


5264.80 


3.230 


-3.23 


7.61 


7.46 


71 


Fe II 


5414.07 


3.221 


-3.48 


7.37 


7.28 


45 


Fe II 


5525.12 


3.267 


-3.94 


7.35 


7.33 


22 



Table 2 

Spectral Lines Selected for the Iron Abundance Analysis of the Sun 



Species 


Wavelength 


EP 


log 5/ 


log e 3D 


loge ID 


W A (3D) 




(A) 


(cV) 


(dex) 


(dex) 


(dex) 


(mA) 


Fe I 


4389.25 


0.052 


-4.55 


7.36 


7.45 


75 


Fe I 


4445.47 


0.087 


-5.38 


7.29 


7.41 


41 


Fe I 


5044.21 


2.851 


-2.02 


7.37 


7.40 


75 


Fe I 


5247.05 


0.087 


-4.98 


7.43 


7.54 


68 


Fe I 


5250.21 


0.121 


-4.90 


7.35 


7.47 


67 


Fe I 


5253.46 


3.283 


-1.57 


7.38 


7.40 


79 


Fe I 


5329.99 


4.076 


-1.22 


7.44 


7.47 


58 


Fe I 


5412.79 


4.434 


-1.72 


7.34 


7.36 


18 


Fe I 


5491.83 


4.186 


-2.19 


7.41 


7.43 


13 


Fe I 


5525.54 


4.230 


-1.08 


7.32 


7.35 


53 


Fe I 


5661.35 


4.284 


-1.76 


7.37 


7.39 


23 


Fe I 


5701.54 


2.559 


-2.14 


7.48 


7.52 


87 


Fe I 


5705.46 


4.301 


-1.36 


7.35 


7.38 


38 


Fe I 


5778.45 


2.588 


-3.44 


7.35 


7.43 


22 


Fe I 


5784.66 


3.396 


-2.53 


7.36 


7.41 


27 


Fe I 


5855.08 


4.607 


-1.48 


7.39 


7.41 


23 


Fe I 


5956.69 


0.859 


-4.50 


7.33 


7.47 


54 


Fe I 


6082.71 


2.223 


-3.55 


7.37 


7.47 


35 


Fe I 


6136.99 


2.198 


-2.93 


7.37 


7.45 


65 


Fe I 


6151.62 


2.176 


-3.37 


7.47 


7.58 


51 


Fe I 


6173.34 


2.223 


-2.88 


7.43 


7.52 


69 


Fe I 


6200.31 


2.608 


-2.37 


7.42 


7.51 


74 


Fe I 


6219.28 


2.198 


-2.45 


7.47 


7.50 


91 


Fe I 


6240.65 


2.223 


-3.17 


7.28 


7.38 


49 


Fe I 


6265.13 


2.176 


-2.54 


7.46 


7.53 


88 


Fe I 


6271.28 


3.332 


-2.70 


7.37 


7.42 


24 


Fe I 


6280.62 


0.859 


-4.39 


7.43 


7.51 


65 


Fe I 


6297.79 


2.223 


-2.64 


7.32 


7.41 


75 


Fe I 


6322.69 


2.588 


-2.47 


7.56 


7.64 


77 


Fe I 


6481.87 


2.279 


-3.01 


7.47 


7.55 


64 


Fe I 


6498.94 


0.958 


-4.69 


7.41 


7.56 


46 


Fe I 


6574.23 


0.990 


-5.02 


7.39 


7.54 


28 


Fe I 


6581.21 


1.485 


-4.68 


7.28 


7.39 


18 


Fe I 


6593.87 


2.433 


-2.37 


7.48 


7.54 


87 


Fe I 


6609.11 


2.559 


-2.66 


7.45 


7.53 


67 


Fe I 


6625.02 


1.011 


-5.34 


7.36 


7.51 


15 


Fe I 


6667.71 


4.584 


-2.11 


7.49 


7.51 


9 


Fe I 


6699.14 


4.593 


-2.10 


7.44 


7.46 


8 


Fe I 


6739.52 


1.557 


-4.79 


7.25 


7.35 


12 


Fe I 


6750.15 


2.424 


-2.59 


7.44 


7.54 


76 


Fe I 


6793.26 


4.076 


-2.33 


7.37 


7.40 


13 


Fe I 


6804.27 


4.584 


-1.81 


7.40 


7.42 


14 


Fel 


6837.01 


4.593 


-1.69 


7.39 


7.42 


17 


Fe I 


6854.82 


4.593 


-1.93 


7.48 


7.50 


13 


Fe I 


6945.21 


2.424 


-2.45 


7.43 


7.51 


83 


Fe I 


6971.93 


3.018 


-3.34 


7.33 


7.40 


13 


Fe I 


6978.85 


2.484 


-2.45 


7.45 


7.49 


80 


Fe I 


7189.15 


3.071 


-2.77 


7.48 


7.54 


40 


Fe I 


7401.69 


4.186 


-1.60 


7.47 


7.51 


42 


Fe I 


7723.21 


2.279 


-3.62 


7.51 


7.62 


41 



Table 2 — Continued 



Species 


Wavelength 


EP 




log e 3D 


loge ID 


W A (3D) 




(A) 


(cV) 


(dex) 


(dex) 


(dex) 


(mA) 


Fe II 


4576.33 


2.844 


-2.91 


7.40 


7.35 


66 


Fe II 


4620.51 


2.828 


-3.19 


7.30 


7.28 


52 


Fe II 


4656.98 


2.891 


-3.58 


7.32 


7.29 


34 


Fe II 


5234.62 


3.221 


-2.23 


7.50 


7.46 


86 


Fe II 


5264.80 


3.230 


-3.23 


7.58 


7.56 


47 


Fe II 


5414.07 


3.221 


-3.48 


7.35 


7.33 


28 


Fe II 


5525.12 


3.267 


-3.94 


7.34 


7.33 


12 


Fe II 


5627.49 


3.387 


-4.10 


7.44 


7.42 


9 


Fe II 


6432.68 


2.891 


-3.51 


7.39 


7.37 


41 


Fe II 


6516.08 


2.891 


-3.38 


7.54 


7.51 


54 


Fe II 


7222.39 


3.889 


-3.37 


7.63 


7.62 


20 


Fe II 


7224.48 


3.889 


-3.29 


7.53 


7.52 


20 


Fe II 


7449.33 


3.889 


-3.07 


7.25 


7.23 


18 


Fe II 


7515.83 


3.903 


-3.45 


7.47 


7.46 


13 


Fe II 


7711.72 


3.903 


-2.45 


7.38 


7.35 


47 
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Table 3 

Radiative Lifetimes for the Different Levels of Fe II 



Term 


J 


Bicmont et al 


Guo et al 


Schnabel et al 
(ns) 


Adopted 


Term average 


6 n 


a /o 

9/2 
7/2 
5/2 
3/2 
1/2 


O T A A A£? 

3.70 0.06 
3.68 0.07 
3.63 0.08 
3.83 0.10 
3.76 0.10 




3.64 0.09 
3.70 0.05 
3.73 0.07 


O 7A A A/^ 

3.70 0.06 
3.66 0.07 
3.68 0.05 
3.76 0.07 
3.76 0.10 


3.71 cr=0.05 


«F 


11/2 

9/2 

7/2 

5/2 

3/2 


3.19 0.04 
3.24 0.06 
3.26 0.10 

3.33 0.09 

3.34 0.10 






3.19 0.04 
3.24 0.06 
3.26 0.10 

3.33 0.09 

3.34 0.10 


3.27 cr=0.06 


6p 


7/2 
5/2 
3/2 


3.73 0.05 
3.83 0.07 


3.73 0.06 
3.79 0.12 
3.71 0.12 




3.73 0.05 
3.81 0.07 
3.71 0.12 


3.75 cr=0.05 


4p 


5/2 
3/2 




3.43 0.09 

3.44 0.11 




3.43 0.09 

3.44 0.11 


3.44 cr=0.01 


4 D 


7/2 
5/2 
3/2 
1/2 




3.02 0.07 
3.10 0.08 


2.97 0.02 

2.90 0.06 

2.91 0.09 


2.97 0.02 
2.97 0.06 
2.91 0.09 
2.95 0.03 a 


2.95 cr=0.03 


4 F 


9/2 
7/2 
5/2 
3/2 




3.87 0.09 
3.63 0.11 
3.75 0.14 


3.55 0.08 


3.87 0.09 
3.63 0.11 
3.60 0.08 
3.70 0.15 a 


3.70 cr=0.15 



a Measurement not available. The term average is adopted. Note that Hannaford et al give 
2.9 ± 0.2 ( 4 D) and 3.7 ± 0.2 ( 4 F) 
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Table 4 

Oscillator Strengths for Fe II Lines in the Optical Region 



Wavelength Term Kroll & Kock 


Whaling Heise & Kock 


Pauls et al 


Adopted 


(A) 


(dex dex ns) 







3002.64 
3163.09 
3170.34 
3183.11 
3185.31 
3186.74 
3192.92 
3193.80 
3196.08 
3210.45 
3213.31 
3227.73 
3255.89 
3277.35 
322.774 
3281.30 
3285.41 
3295.81 
3302.86 
3303.47 
4173.47 
4178.87 
4233.17 
4303.17 
4351.76 
4385.39 
4491.40 
4508.29 
4515.33 
4522.62 
4549.46 
4555.88 
4576.33 
4582.83 
4583.83 
4620.51 
4629.34 
4656.97 
4923.92 
5169.03 
5197.57 
5234.66 
5262.47 
5264.79 



l P: 



5/2 
7 5/2 
^1/2 
7 5/2 
7 3/2 
^3/2 
^5/2 
^1/2 
r 7/2 
3/2 
^5/2 
D 7 / 2 

6 D 7/2 
6 D 9/2 
i D 7 / 2 
6 D 5/2 
6 D 1/2 
6 D 3/2 

6 D7/2 

6 D 1/2 
4 D 5/ 2 
F7/2 
D7/2 

3/2 
5/2 
1/2 
3/2 



^1/2 

7 5/2 
^3/2 



F; 



5/2 



4p 



4 D 

6R 



7/2 
5/2 
7/2 
7/2 

^7/2 
^9/2 
5/2 
'3/2 
7/2 
3/2 
5/2 



5/2 
3/2 



-0 


.93 


0.07 


3.9 








-0.87 0.07 


-2. 


,82 


0.07 


4.0 




-2.82 0.07 4.0 a 


-2.49 0.08 3.4 


-2.78 0.07 
-2.42 0.08 


-2. 


.10 


0.07 


4.0 


-2.04 0.11 4.0 






-2.04 0.06 


-2. 


.78 


0.07 


4.1 








-2.74 0.07 


-1 


.67 


0.04 


3.2 


-1.71 0.11 3.4 






-1.63 0.04 


-1 


.95 


0.11 


3.2 






-1.92 0.07 3.4 
-1.75 0.07 3.4 


-1.88 0.06 
-1.69 0.07 


-1 


.73 


0.11 


3.9 


-1.66 0.11 3.9 






-1.67 0.08 


-1 


.69 


0.04 


3.2 


-1.79 0.11 3.4 


-1.69 0.04 3.2 a 




-1.66 0.04 


-1 


.27 0.07 


3.2 






-1.31 0.07 3.4 


-1.25 0.05 


-1 


Ofi 
. uu 


1 1 


3.7 


-1.13 0.04 3.7 






-1.02 0.04 


-2. 


.52 


0.04 


4.2 


-2.50 0.22 4.0 


-2.52 0.04 4.2 a 




-2.46 0.04 


-2. 


'-id 


04 


3.9 


-2.47 0.11 3.9 
-1.13 0.07 3.7 






-2.31 0.04 
-1.03 0.07 


-2 


.69 


0.04 


4.1 




-2.69 0.04 4.1 a 




-2.65 0.04 


-2 


.87 


0.04 


3.9 


-2.12 0.06 4.0 






-2.82 0.04 


-2 


9(1 


04 


4.0 




-2.90 0.04 4.0 a 




-2.87 0.04 


-3 


.51 


0.07 


4.2 








-3.45 0.07 


-2 


.70 


0.11 


3.9 








-2.68 0.11 


-2 


.18 


0.07 


3.2 






-2.77 0.14 3.4 


-2.50 0.08 


-2 


.48 


0.11 


3.9 








-2.45 0.11 


-1 


.91 


0.07 


3.7 


-2.00 0.11 3.7 


-1.91 0.07 3.7 a 




-1.84 0.06 


-2 


.49 


0.07 


3.2 




-2.65 0.07 3.2 




-2.61 0.07 


-2 


.10 


0.07 


3.2 




-2.10 0.07 3.2 a 


-1.99 0.08 3.4 
-2.75 0.06 3.4 


-2.03 0.05 
-2.69 0.06 


-2 


.70 


0.11 


4.1 






-2.58 0.09 3.4 


-2.66 0.11 
-2.52 0.09 


-2 


.41 


0.11 


4.0 








-2.36 0.11 


-2 


.03 


0.11 


3.2 








-1.99 0.11 


-1 


.75 


0.07 


3.2 






-2.33 0.30 3.4 


-1.87 0.08 


-2 


.29 


0.11 


3.9 




-2.94 0.09 3.2 


-2.97 0.16 3.4 


-2.26 0.11 
-2.91 0.07 


-3 


.10 


0.11 


3.9 








-3.07 0.11 


-1 


.84 


0.07 


3.7 


-2.02 0.22 3.7 


-3.29 0.08 3.7 




-1.77 0.07 
-3.19 0.07 


-2 


.33 


0.11 


4.3 




-3.61 0.10 3.2 




-2.28 0.11 
-3.58 0.08 


-1 


.24 


0.11 


4.0 


-1.32 0.11 4.0 






-1.26 0.08 


-0 


.87 


0.11 


3.8 








-0.86 0.11 


-2 


.10 


0.11 


4.1 








-2.06 0.11 


-2 


.05 


0.11 


4.0 




-2.27 0.10 4.0 d 
-3.06 0.10 3.2 d 
-3.27 0.08 3.2 




-2.23 0.08 
-3.03 0.10 
-3.23 0.05 
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Table 4 — Continued 



Wavelength Term Kroll & Kock Whaling 


Hcise & Kock 


Pauls et al 


Adopted 


(A) 


(dex dex ns) 







5276. 


.00 


4 F 7 / 2 


5316. 


.62 


F9/2 


5316. 


.78 


4 D 5/2 


5325. 


.55 


A F 7/2 


5414. 


,08 


-D7/2 


5525. 


.13 




5607. 


,12 


& D 7/2 


5627. 


.49 


4 F 5/2 


6369 


,42 


6 D 3/2 


6432. 


,67 


6 D 5/2 


6516. 


,07 


6 D 7/2 


7222. 


.39 


4 D 1/2 


7224. 


,49 


A D l/2 


7449. 


,34 


4 As/2 


7515. 


,83 


Z 5/2 


7711. 


,73 


^7/2 



-1.94 0.11 3.9 
-1.85 0.11 4.3 

-2.60 0.11 3.9 

















-1.91 





.11 
















-1.85 


0. 


.11 














-2.80 0.12 3.4 


-2.74 





.12 
















-2.57 





,ll b 


-3 


,58 


0, 


.09 


3. 


.7 




-3.48 





.08 


-4 


,04 





.11 


3 


.7 




-3.94 





.09 


-3 


,83 





.07 


4. 


.2 




-3.77 





.07 


-4 


.14 


0. 


.09 


4. 


.0 




-4.10 





.07 


-3 


.50 


0. 


.12 


4. 


.0 




-3.47 





.12 


-3 


.55 


0. 


.08 


4. 


.1 




-3.51 





.06 


-3 


.44 


0, 


.07 


4. 


.2 




-3.38 





.05 














-3.43 0.07 3.4 


-3.37 





.07 














-3.35 0.06 3.4 


-3.29 





.06 


-3 


.10 


0, 


.11 


3. 


.2 




-3.07 





.09 


-3 


.41 


0. 


.12 


3 


.2 


-3.53 0.08 3.4 


-3.45 





.07 


-2 


.55 


0, 


.08 


3. 


.7 




-2.45 





.07 



a HK lines taken as reference from KK 

b This line has log gf — —3.2 as computed by Kurucz. Only with a value similar to that we can 
reconcile the iron abundance from this line with that from the other Fe II lines 

c Comparison with the solar spectrum shows that a value closer to -3.38 is to be preferred 
d Line classified as a possible blend or a misidcntification by Heise & Kock 



